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Easy access to next generation sequencing has enabled the rapid analysis of complex microbial populations. To take full 
advantage of these technologies, animal models enabling the manipulation of human microbiomes and the study of 
the impact of such perturbations on the host are needed. To this aim we are developing experimentally tractable and 
clinically relevant pig models of the human adult and infant gastro-intestinal tract. The intestine of germ-free piglets was 
populated with human adult or infant fecal microbial populations, and the piglets were maintained on solid or milk diet, 
respectively. Amplicons of 16S rRNA V6 region were deep-sequenced to monitor to what extent the transplanted human 
microbiomes changed in the pig. Within 24 hours of transfer of human fecal microbiome to pigs, bacterial microbiomes 
rich in Proteobacteria emerged. These populations evolved toward a more diverse composition rich in Bacteroidetes and 
Firmicutes. In the experiment where infant microbiome was used, the phylogenetic composition of the transplanted 
bacterial population converged toward that of the human inoculum. A majority of sequences belonged to a relatively 
small number of operational taxonomic units, whereas at the other end of the abundance spectrum, a large number 
of rare and transient OTUs were detected. Analysis of fecal and colonic microbiomes originating from the same animal 
indicate that feces closely replicate the colonic microbiome. We conclude that the pig intestine can be colonized with 
human fecal microbiomes to generate a realistic model of the human Gl tract. 



Introduction 

Intestinal microbiomes play an important role in human health. 
Dysbiotic intestinal bacterial populations have been linked to a 
variety of health conditions, 1 ' 2 even with conditions not directly 
associated with the gastro-intestinal (GI) tract such as asthma, 3 
diabetes, 4 disorders of the immune system, 5,6 arthritis 7 and preg- 
nancy. 8 Research on intestinal microbiota and their interaction 
with the host often focuses on human subjects or uses germ- 
free mice inoculated with human fecal microbiota. The search 
for associations between human phenotypes and microbiome is 
hampered by the many confounding variables inherent to human 
cohorts and by the small number of subjects typically included 
in such studies. Humans are genetically diverse, are exposed to 
many environmental factors and consume different diets, all fac- 
tors which directly or indirectly impact the intestinal ecology. To 
circumvent these limitations and to study the effect of specific 
alterations of the intestinal microbiome, germ-free rodents have 
been populated with microorganisms originating from human 
feces.'" 14 Rodent models, and especially murine models, offer 
many advantages over other species. The fact that many strains 
are inbred and the availability of genetically modified lines 
facilitate research aiming at elucidating the interaction between 
microbiome, genetic background and disease. 15 On the other 
hand, rodent models are limited by the many physiological dif- 
ferences between rodents and primates. Most importantly, mice 
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and rats typically do not display clinical manifestations seen in 
human enteric diseases. 

The need for a clinically relevant model of the human GI 
tract is driving the search for non-rodent models. The advantage 
of pigs as an alternative model of human diseases lies in similar 
clinical manifestations and their susceptibility to many enteric 
pathogens afflicting humans. Efforts to develop porcine models 
benefit from extensive research on the use of pigs as source of tis- 
sues and organs for xenotransplantation. 16 This research expands 
our knowledge on this species and has led to the production of 
transgenic pigs. 17 It has also promoted the production of immune 
reagents to study the porcine immune response. 

Piglets derived by C-section are colostrum-deprived and 
hence, unlike humans and rodents, are born agammaglobulin- 
emic, with no interfering maternal antibodies. As with human 
infants, the pig's immune system at birth is fully developed but 
not fully functional, maturing within 7-10 d. Relevant features 
of the immune system include a distribution of Peyer's patches 
which is similar to that seen in humans and play a central role 
in antigen sampling by M cells. Pigs are also immunologically 
similar to humans in their pattern of lymphocyte distribution in 
the body and at the mucosal surfaces. Intra-epithelial lympho- 
cytes are also broadly similar to those of humans, mice and rats, 
with a majority of cytotoxic suppressor T cells and fewer T helper 
cells. The presence of gut microflora in these animals acceler- 
ates gut maturity and promotes a robust immune response. 18 " 20 
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Figure 1. Phylum-level taxonomy of fecal and colonic bacterial populations from pigs inoculated with human fecal microbiomes. Experiments 1-3 
refer to the three experiments described in the text and in Table 1. Phyla are color-coded as shown in the key. A total of 10" V6 16S sequence reads per 
sample were classified. The taxonomic time series of two pigs from experiments 1 and 3 are shown. A single time series is shown for experiment 2. The 
x axis shows pig age in days. Bars labeled "H" represent human inocula; samples obtained directly from the colon are labeled "gut." The two left-most 
bars in experiment 2 are from independent replicate analyses of the infant fecal sample used as inoculum. 



Experiments have been conducted in the pig to study the effect 
of gut microflora on the immune system, both in conventional 
and germ-free animals. 21 The protective effect of probiotics on 
subsequent colonization with pathogenic bacteria has been dem- 
onstrated in this species. 20,22 

From a practical point of view, piglets learn to drink sterile 
milk formula directly from a trough. At three weeks of age they 
can be weaned onto solid diet. Litters of 8-12 piglets provide 
a good sample size. Newborn piglets weigh about 1 kg, which 
facilitates handling, bleeding, surgery, repeated sampling, mea- 
suring of temperature, heartbeat, respiratory rates, collection 
of intestinal content for microbiome analyses and of intestinal 
epithelial cells for analyzing the transcriptional response. Piglets 
derived prematurely can be used to study complications of pre- 
term births. 

Here we describe results obtained from three experiments 
comprising eight germ-free piglets inoculated with human adult 
or infant fecal microbiome. The transplanted microbial popula- 
tions were monitored for up to 35 d. Evidence supporting the 
successful transplantation of human microbiota is presented. The 
implications of these observations for the development of better 
models of the human GI tract are discussed. 

Results 

Overview. A summary of the pig experiments is shown in Table 1 . 
A total of three experiments were performed; two experiments 



using fecal extracts from human adults and one with fecal extract 
from a 3 mo old breast-fed baby. Small portions of these extracts 
were stored at -80°C to enable multiple experiments with the 
same stock. In experiment 1, adult microbiota was transplanted 
into two pigs 8 d of age maintained on Similac milk replacer for 
the 20-d duration of the experiment. Experiment 2 and 3 were 
designed to emulate the infant and adult GI tract, respectively. 
In experiment 2, infant microbiota was inoculated on day 5 or 
30 into four germ-free piglets maintained on Similac, whereas in 
experiment 3 two pigs were inoculated with adult microbiota on 
day 23 of age, three days after weaning them onto solid porcine 
diet. 

Phylogeny. The phylum-level classifications of time series 
originating from five pigs from experiments 1—3 are shown in 
Figure 1. In general, a reduction in the proportion of Proteobacteria 
and expansion of Firmicutes was observed. This trend was partic- 
ularly apparent in experiment 2, which likely reflects the relative 
high proportion of Proteobacteria in the infant inoculum, and the 
fact that fecal samples collected 24 h post-inoculation were ana- 
lyzed in this experiment. The phylogenetic profile of the earliest 
sample from three other pigs included in this experiment were 
also rich in Proteobacteria (98%, 78%, 84%). The histograms 
for those animals are not shown because of the short time series 
(see Table 1). Samples collected on the day following inoculation 
were not available for experiment 1 and 3. We noticed in experi- 
ments 2 and 3 a close resemblance in the taxonomic profile of 
the fecal and colonic microbiome, labeled "gut" in Figure 1. The 
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Table 1. Summary of pig experiments 





experiment 1 




experiment 2 




experiment 3 


inoculum 


adult 




infant 




adult 


pig ID 


4 5 
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1 2 


age inoculated 3 


8 8 
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5 30 


30 


23 23 


diet 


Similac Similac 


Similac 


Similac Similac 


Similac 


solid solid 


day euthanized 


28 28 


41 


15 41 


41 


36 36 



a Age in days. 
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Figure 2. Rank-abundance plots of 16S sequence reads. Experiments 1-3 refer to the three experiments described in the text. Green indicates human 
sample. Samples collected from a same animal on subsequent days are coded with the same color as indicated in the keys. Some loss of diversity is 
apparent, particularly in experiment 1. 



minor differences between fecal and colonic profile are of the 
same magnitude as between replicate human profiles shown in 
experiment 2 and are thus likely to represent technical variability. 
Unclassified sequence reads were particularly abundant in experi- 
ment 3. As further discussed below, this observation is probably 
related to sequencing strategy based on 60-nt reads. Providing a 
more detailed view of the evolution of the transplanted human 
fecal microbiomes, genus-level taxonomic classifications for the 
same five pigs as shown in Figure 1 are shown in Table SI. The 
table shows raw counts of sequences assigned to each genus. 

Analysis of sequence evenness and diversity. To compare 
the evenness of intestinal microbiota before and after trans- 
plant into pigs, rank-abundance plots for each experiment were 
drawn. These plots showed some loss of evenness, particularly 
in experiments 1 and 3, where fecal microbiome from adults 
was transferred to milk-fed pigs (Fig. 2). The pig microbiomes 
which showed little or no loss of evenness were in experiment 2, 
in which the microbiome of a breast-fed child was transplanted to 
pigs fed milk replacer. In experiment 3, adult microbiome trans- 
planted into weaned pigs experienced some loss of evenness. 

Alpha diversity was estimated using the Shannon index 
(Fig. 3). Consistent with the rank-abundance plots, experiment 2 
shows that the pig microbiomes recovered to their original diver- 
sity after about 5 weeks in the pig. In contrast, in experiments 
1 and 3 the pig microbiomes lost diversity, a trend that did not 
reverse itself over the duration of the experiment. As diversity in 
experiment 2 increased toward the end of the 4l-d time series, it 
is conceivable that the same would have occurred in experiments 



1 and 3 had the time series been extended to 41 d. The change 
in diet experienced by the microbiome in these two experiments 
could also have contributed to a loss in bacterial diversity. 

Analysis of UniFrac distances. PCoA was used to visualize 
weighted pairwise UniFrac distances. A separate analysis was per- 
formed for each of the 3 experiments to assess the divergence of the 
microbiome following transplant into pigs. To avoid data point 
compression, data from samples collected within 24 h of inocula- 
tion were excluded, as these populations were characterized by 
a high abundance of Proteobacteria. To visualize the trajectory 
of the transplanted fecal bacterial populations, data points from 
subsequent days were linked (Fig. 4). The PCoA plots illustrate 
the extent to which the transplanted fecal microbiomes evolved 
in the pig GI tract over a period of up to 35 d. Consistent with 
the evenness and diversity analyses shown in Figures 2 and 3, the 
data points from experiment 1 reveal a relatively large distance 
between inoculum and transplanted populations. Pig microbi- 
omes collected at four time points over a period of 15 d did not 
converge toward the composition of the human inoculum. These 
experiment 1 trajectories contrast with those observed in experi- 
ment 2. In this experiment, bacterial populations transplanted 
into four pigs converged toward the human inoculum, consis- 
tent with a similar diversity of the human inoculum and the pig 
microbiome 41 d post-inoculation. The trajectories for pigs 8 and 
10 are only three and two days long, respectively, which makes it 
difficult to assess whether the microbiome would ultimately have 
converged toward the human inoculum. Pig 8 was euthanized 
on day 15 due to deteriorating health and for pig 10 only two 
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Figure 3. Diversity of pig intestinal microbiome OTUs. Diversity was estimated using the Shannon diversity index calculated with the natural log. 
Experiments 1 -3 refer to the three experiments described in the text. The x axis shows pig age in days. Black bars indicate human inocula. Bars are 
shaded according to pig; experiments 1 and 3 included two pigs each, shown with white and gray bars, respectively. Four pigs were used in experi- 
ment 2, as indicated with white (pig 7), light gray (pig 8), dark gray (pig 9) and hatched bars (pig 10). 



fecal samples and the gut content were sequenced. In experiment 
3, adult microbiome transplanted into 23-d old pigs fed a solid 
diet did not follow an obvious trajectory in the course of the time 
series lasting 12 d. Weighted UniFrac distances from the inocu- 
lum 2 d post-inoculation were 0.476 and 0.408 for pig 1 and pig 2, 
respectively. On day 35 of age weighted distances from the inocu- 
lum were almost unchanged at 0.421 and 0.451, for pig 1 and 2, 
respectively. 

A significant advantage of animal models over human subjects 
is the possibility to easily collect samples directly from the colon. 
We took advantage of this possibility to compare microbiomes 
originating from the feces excreted prior to euthanasia and col- 
lected from the colon a short time after. Samples from the colon 
were obtained from four of the eight animals examined here, as 
indicated with an open circle in Figure 4 and as labeled with "gut" 
in Figures 1 and 3. In all cases the colonic samples closely resem- 
bled the fecal microbiomes from the same animal in terms phy- 
logeny (Fig. 1), diversity (Fig. 3) and UniFrac distance (Fig. 4). 

Discussion 

An ideal model of the human GI tract faithfully replicates the 
ecology of the gut. Such a model would be susceptible to dis- 
eases afflicting humans, develop similar symptoms, mimic the 
adults or infant GI tract and would be experimentally tractable 
and affordable. The pig model we present above fulfills several 
criteria; specifically it is experimentally tractable and it is a clini- 
cally relevant model of the human GI tract. We aim to use this 
model to explore the interaction between microbiome and host in 
healthy and diseased animals. Such a model will enable research 
on how to cure or alleviate disease of the GI tract. 

For obvious reasons, much of the research community is using 
rodents as a model to study the human GI tract. 11,13,23 " 25 Although 
the pig model we describe here is not perfect, we believe that it 
offers significant advantages over rodents for basic and transla- 
tional research. Foremost is the physiological similarity between 
the human and porcine GI tract. Pigs have not only been shown 
to be susceptible to many enteric human pathogens, 26 " 32 but pres- 
ent symptoms characteristic of human infections. 33 The results 



described here showing that a human intestinal ecosystem can be 
transplanted into pigs expands the potential role of pigs in research 
beyond the germ-free or gnotobiotic model. The data presented 
here indicate that pigs are a particularly good model for study- 
ing the development of the infant microbiome. For instance, the 
infant model would enable studies on how the exposure to vari- 
ous components of the human microbiome affects development 
of the innate immune system. The research we have presented 
paves the way for experiments aiming at elucidating a more com- 
plex and realistic model of the GI tract involving the host, the 
microbiome and enteric pathogens. We intend to perform future 
research with the same human donor samples which have been 
preserved at -80°C. Although we have not tested to what extent 
storage conditions have affected the composition of the microbi- 
ome, we believe that the advantage of using a standardized stocks 
of adult and infant microbiomes overrides any concerns about the 
impact of storage on the microbial population. 

Our sequencing strategy is based on the Illumina platform and 
differs from more commonly used pyrosequencing methods. 34 
The main difference to pyrosequencing is the shorter length of 
sequence reads. The Illumina platform was chosen for practical 
reasons, specifically because its higher throughput enables the 
analysis of an essentially unlimited number of samples collected 
from multiple animals over an entire experimental period. 
The fact that the V6 amplicon sequences are compatible with 
phylogenetic classification and diversity analysis is due to several 
factors. In the first place, the V6 loop targeted in our protocol has 
a higher information entropy than any other region in the 16S 
gene. 35 Second, we used a 16S database of reference sequences 
trimmed to 100 nt centered on the V6 loop. This approach 
improves classification. A systematic evaluation of the effect of the 
sequence database showed that customized reference databases, 
as used here, can reduce the number of unclassified reads. 36 We 
verified the observation by Werner et al. using the conventional 
full-length RDP template database 37 or the V6 database. In 
agreement with Werner et al. 36 we found a reduction in the 
proportion of unclassified reads when using a trimmed template 
file. Experiments in which classification based on the V6-V9 
region was compared with a truncated template corresponding to 
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Figure 4. Principal coordinate analysis of 16S sequences from human fecal microbiota before and after transplant into germ-free pigs. The analysis is 
based on pairwise weighted UniFrac distances. Human microbiota are shown in green. Data points from each pig, color-coded as shown in the keys, 
are connected in chronological order. Numbers indicate age in days of first and last sample. Full symbols indicate microbiota from the feces, empty 
symbols samples recovered from the colon. The proximity of the gut microbiota and the fecal microbiota obtained on the final day of the experiment 
(day 41, 38 and 35, respectively) indicate a close similarity between colon and fecal microbiota. 



the V6 region support the view that a short hypervariable region 
can be used instead of longer regions. This analysis revealed that 
classification was more likely to be impacted by primer bias than 
by read length. 38 

A system-wide analysis of the gut environment and its 
response to perturbation will require expanding the analyses 
beyond 16S amplicon sequencing to the microbial metagenome 
and the transcriptome of the gut epithelium. Including the epi- 
thelial transcriptome will reveal if and how the host responds to 
perturbations of the intestinal ecosystem, for instance as a result 
of a change in diet, as a result of infection or treatment with anti- 
microbials. Surprisingly, we found only one report of such a study 
which used microarrays to compare the epithelial transcriptome 



of germ-free piglets and piglets populated with porcine micro- 
biome. This study observed transcriptional changes of genes 
involved in epithelial cell turnover, mucus synthesis and in the 
regulation of the immune system. 39 Because of the larger size of 
the pig intestine, as compared with the mouse, we anticipate that 
the pig model will enable a more comprehensive analysis of the 
gut environment, including microorganisms which live in close 
association with the epithelium 40 ' 41 and may not be adequately 
represented in fecal samples. 

Although the number of experiments we have performed to 
date is small and does not support robust inferences on the puta- 
tive impact of diet on microbiome diversity, our observations (see 
Figs. 2 and 3) are consistent with surveys in humans 42,43 and 
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experiments in rodents. 11,44 ' 45 If confirmed with additional experi- 
ments in pigs, the putative impact of diet on the composition of 
the microbiome would support the view that perturbations of the 
human GI tract environment can be replicated in microbiomes 
transplanted into pigs. 

In contrast to a study of the evolution of native microbiome in 
mice over a period of one year, 25 our study did not identify a mea- 
surable stabilization on the microbiome over time. In our experi- 
ments, weighted Unifrac distances typically decreased during the 
first few days following inoculation, but no obvious trend emerged 
thereafter. We hypothesize that the relatively short duration of the 
pig experiments was not sufficient for the gut microbiome to stabi- 
lize. Because our animals are housed in microbiological isolators 46 
for the duration of the experiment, extending the experiment 
beyond 5 or 6 weeks is not feasible. Experiments are needed to 
assess whether transferring pigs colonized with human microbi- 
ome to a normal environment will affect the composition of the 
microbiome. If, once established in the pigs' GI tract, the human- 
ized microbiome becomes resistant to invasion by environmental 
bacteria, experiments of longer duration would become possible. 
Such experiments could enable research on slow-evolving or 
chronic conditions and reveal to what extent intestinal dysbiosis is 
associated with chronic conditions. 

Materials and Methods 

Animals and human fecal samples. Piglets were derived by 
C-section and housed in sterile isolators. 46,47 Piglets were fed three 
times daily with human infant milk formula (Similac, Abbott). 
They were checked for bacterial contamination with daily aero- 
bic culture on Brucella Blood Agar. In one experiment (experi- 
ment 3) pigs were weaned on day 20 by feeding ad libidum sterile 
Laboratory Porcine Grower Diet with 16% protein (LabDiet). 

Feces were collected from 10 adult human donors (5 males 
and 5 females) aged 50 to 70 y and from a 3-mo old breast-fed 
healthy baby. To minimize the loss of anaerobic microorganisms, 
feces were transferred into completely filled airtight containers. 
Samples were homogenized in 9 volumes of reduced PBS con- 
taining 10% glycerol. The homogenates were filtered through a 
425-|xm pore mesh, dispensed into 15-ml centrifuge tubes and 
stored at -80°C. The samples from the adult donors were com- 
bined into a single sample. 

On the day of inoculation with human fecal microbiota 
(Table 1), piglets were first given 10 ml of 0.2 M carbonate buffer 
pH 9.5 orally, followed by 3 ml of stool homogenate with a feeding 
needle inserted into the esophagus. The animals were monitored 
for the duration of the experiment for signs of disease, including 
diarrhea, dehydration, dyspnea, weakness, lethargy, anorexia, etc. 
Weights were recorded every other day and fecal samples collected 
for bacterial culture, for measuring cytokines and for extracting 
DNA for high-throughput sequencing. Piglets were euthanized 
with an intramuscular injection of 1 ml 100 mg/ml ketamine and 
0.1 ml 100 mg/ml xylazine per kg body weight, followed by an 
intra-cardial injection of 1 ml of 390 mg/ml SomnaSol/kg (Butler 
and Schein). Following euthanasia, gut contents were recovered 
from the spiral colon by gently squeezing the colon. 



DNA extraction and 16S library preparation. DNA was 

extracted from 200 ixl of fecal slurry. Fecal samples were first sub- 
jected to two cycles of freeze-thawing (-80 C C/37°C). DNA was 
extracted from these samples with the HighPure PCR Template 
Preparation kit (Roche Diagnostics). The final DNA extract was 
dissolved in 20 (Jtl buffer. A volume of 1 |xl of DNA was ampli- 
fied in a primary PCR with primers flanking the approximately 
60-nucleotide variable region of the V6 loop of the bacterial 16S 
rRNA gene. 48 The following primers were used: 5-CAA CGC 
GAA GAA CCT TAC C-3' and 5-CGA CAG CCA TCG ANC 
ACC T-3'. The primary amplification consisted of 15 cycles of 
94°C for 30 sec, 55°C of 30 sec and 68°C for 90 sec. A sec- 
ondary PCR was used to amplify 1 ixl of primary amplicon and 
generate a library of PCR products fused to Illumina adaptors 
and tagged with a 6-nucleotide barcode unique to each sample. 
The downstream (reverse) primer used in the secondary reaction 
incorporates in 5' to 3' orientation the Illumina flow-cell binding 
sequence, the barcode, the complement of the standard Illumina 
multiplex index read primer and the downstream conserved 
sequence flanking the V6 region. The secondary PCR consisted 
of three cycles of 94°C for 30 sec, 55°C of 30 sec, 68°C for 90 sec 
followed by 11 cycles of 94°C for 30 sec, 68°C for 30 sec. ExoSap 
was used to remove excess primers from the final PCR and the 
quality of the amplicons verified by electrophoresis on a 2% aga- 
rose gel. The expected amplicon size is approximately 200 nt. 
The pooled barcoded libraries were sequenced in a HiSeq2000 
Illumina sequencer at the Tufts Genomics core facility (tucf.org). 
An average of 3.96 million 100-nucleotide (nt) reads were obtained 
for each sample (range 1.4 x 10 6 — 8.1 x 10 6 reads) (Table 2). 
Quality control for Illumina sequencing included a PhiX phage 
library sequenced with each reaction. This control was used to 
estimate phasing and matrix correction parameters for base call- 
ing. It also provided a reliable statistic of run error rate. 

Data analysis. For most analyses the open-source software 
mothur version 1.25 was used. 49 Sequences were initially parsed 
according to barcode and trimmed to eliminate the downstream 
primer sequence, leaving sequences of -60 nt length in average 
(2.5 percentile = 56 nt; 97.5% percentile = 63 nt). A random sub- 
sample of 10,000 sequences was used for downstream analyses. 
Sequences were classified using the RDP classifier program imple- 
mented in mothur using the bayesian method. 50 To optimize the 
taxonomic classification, 36 a training set template file trimmed 
to a 100-nt region centered on the V6 region was generated from 
the 8422-sequence RDP training set version 6. 37 The minimum 
bootstrap value for taxonomic assignment was set at 70%. 

The phylogenetic distance between samples was quantified 
with the weighted UniFrac distance metric 51 as implemented in 
mothur. Random subsamples of 10,000 sequences used for this 
analysis were screened to identify sequences which did not meet 
the following quality criteria: minimum length 50 nt, maximum 
length 70 nt, start and end of sequence within start and end posi- 
tion of 90% of the sequences. Between 92.0 and 98.1% of the 
sequences met these criteria and were included in the compu- 
tation of UniFrac distance. Lower-triangle matrices of pairwise 
UniFrac distances were imported into GenAlex 52 and distances 
between samples visualized using the Principal Coordinate 
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Table 2. Summary of sequencing results 



samples 




adult (1,3)» 


infant b (2) 


pig (i) 


pig (2) 


pig (3) 


number of samples sequenced 


1 


1 


8 


23 


22 


number of pigs 






2 


4 


2 


sequences meeting quality criteria 


9713 


9200 


9814(81.4) c 


9306 


9581 (359) 


unique sequences 


16.7% 


8.4% 


9.1% (1.2) 


5.8% (1.4) 


13.2% (0.9) 



Experiment number is indicated in parenthesis (see text); pooled adult fecal sample was used for experiments 1 and 3; b mean of two replicates (two 
samples amplified and sequenced independently); standard deviation is shown in parenthesis. 



Analysis (PCoA) calculator. The evenness of human and pig 
samples were compared using mothur's rarefaction. single com- 
mand and visualized on rank-abundance plots. 53 The Shannon 
diversity index was used to estimate a diversity. 
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